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ABSTRACT 

We present ground-based limb polarization measurements of Jupiter and Saturn consisting 
of full disk imaging polarimetry for the wavelength 7300 A and spatially resolved (long slit) 
spectropolarimetry covering the wavelength range 5200 to 9350 A. 

For the polar region of Jupiter we find for A = 6000 A a very strong radial (perpendicular 
to the limb) fractional polarization with a seeing corrected maximum of about +11.5 % in 
the South and +10.0 % in the North. This indicates that the polarizing haze layer is 
thicker at the South pole. The polar haze layers extend down to 58° in latitude. The 
derived polarization values are much higher than reported in previous studies because of the 
better spatial resolution of our data and an appropriate consideration of the atmospheric 
seeing. Model calculations demonstrate that the high limb polarization can be explained 
by strongly polarizing {p ^ 1.0), high albedo {u ~ 0.98) haze particles with a scattering 
asymmetry pararneter o f g ~ 0.6 as expected for aggregate particles of the type described by 



West and Smith! (Il99l[ ). The deduced particle parameters are distinctively different when 



compared to lower latitude regions. 

The spectropolarimetry of Jupiter shows a decrease in the polar limb polarization to- 
wards longer wavelengths and a significantly enhanced polarization in strong methane bands 
when compared to the adjacent continuum. This is a natural outcome for a highly polariz- 
ing haze layer above an atmosphere where multiple scatterings are suppressed in absorption 
bands. For lower latitudes the fractional polarization is small, negative, and it depends only 
little on wavelength except for the strong CH4-band at 8870 A. 

The South pole of Saturn shows a lower polarization {p ^ 1.0 — 1.5 %) than the poles 
of Jupiter. The spectropolarimetric signal for Saturn decrease rapidly with wavelength and 
shows no significant enhancements in the fractional polarization in the absorption bands. 
These properties can be explained by a vertically extended strat ospheric haze region com- 



posed of small particles < 100 nm as suggested previously by iKarkoschka and Tomasko 



(]2005[ ). 



In addition we find in the V- and R-band a previously not observed strong polarization 
feature {p = 1.5 — 2.0 %) near the equator of Saturn. The origin of this polarization signal 
is unclear but it could be related to a seasonal effect. 

Finally we discuss the potential of ground-based limb polarization measurements for the 
investigation of the scattering particles in the atmospheres of Jupiter and Saturn. 
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Introduction 



Reflected light from planets is polarized. Polarimetry therefore provides information on 
the nature and distribution of the scattering pa rticles in the atmospheres of planets comple- 



mentary to other techniques fse e ICoffeenlll979f). As the scattering polarization from extra- 



solar planets can be quite high (IKattawar and Adamslll97ll : ISeager et al.ll2000t IStam et al. 
2OO4J : iBuenzli and Schmidll2009l ). polari metry is also used for the search of extra- solar planets 
with existing and future instruments ( iHough et al.l l2006l : ISchmid et al.l l2006bl ) The polari- 
metric investigations of extra-solar planets also generate renewed interest in the polarization 
properties of solar system planets. These can be used to predict the expected polarization 
signal for extra-solar planets, and thus help in the interpretation of future detections. 

Solar system planets have been frequently observed polarimetri cally until 1990 with 



instruments using single channel (aperture) detectors (e.g. lLeroyll2000l ). However, almost no 
data were taken with "modern" , ground-based imaging polarimeters and spectropolarimeters 
using array detectors. Therefore the polarimetric properties of solar system planets are still 
not well characterized. 

We have therefore started a p rogram of "modern" , grou nd-based polarim e tric ob ser- 
vations of solar system planets. In ISchmid et al.l (l2006al ) and IJoos and SchmidI (l2007al ) we 
described the data for Uranus and Neptune, for which we detected a strong limb polariza- 
tion. In this paper we present imaging polarimetry for Jupiter and Saturn, taken with the 
Zurich imaging polarimeter (ZIMPOL) at the McMath-Pierce solar telescope, and long slit 
spectropolarimetry taken with the EF0SC2 instrument attached to the ESO 3.6m telescope. 

For the outer planets the possible phase angles for ground based observations are very 
limited, and the disk integrated polarization is close to zero due to the back-scattering sit- 
uation. But with spatially resolved observations one can use the limb polarization effect to 
constrain polarimetric properties of the atmosphere. The limb polarization is a well-known 
second order effect occur ring in reflecting at mospheres where Rayleigh-type scattering pro- 
cesses are dominant (e.g. Ivan de Hulstlll980l ). To understand this effect, one has to consider 
a back-scattering situation at the limb of a sphere, where we have locally a configuration of 
grazing incidence and grazing emergence for the incoming and the back-scattered photons, 
respectively. Photons reflected after one scattering are unpolarized, because the scattering 
angle is 180°. Photons undergoing two scatterings travel after the first scattering predom- 
inantly parallel to the surface before being reflected towards us by the second scattering 
process. Photons going up will mostly escape without a second scattering, and photons 
going down have a low probability of being reflected towards us after the second scattering 
but a high probability to be absorbed or to undergo multiple scatterings. Because the po- 
larization angle induced in a single dipole-type scattering process, like Rayleigh scattering. 
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is perpendicular to the propagation direction of the incoming photon (which is often parallel 
to the limb), a net polarization perpendicular to the limb is produced. 

P olarim etric data of Jupiter and Saturn were first published in the pioneering paper of 

LyotI ( 1l929l ). who detected for Jupiter a strong positive polarization of p ~ 5 — 8 % at the 



poles with an orientation perpendicular to the limb. In the disk center he measured a phase 
angle dependent polarization which is essentially zero near opposition and slightly negative 
(parallel to the scattering plane), p ~ —0.4 %, for phase angles around a = 10°. These 
measurements we re confirmed and improved in many later observations using single aperture 



polarimeters (e.g. lDollfuslll957l : iGehrels et al 



and some imaging polarimetry by ICarlson and Lutzl ( 119891 ) . 



1969 



Morozhenkolll973l : iHall and Rilevlll976r ) 



Important results on the polarization of Jupiter were achieved with the Pioneer 10 and 
11 spacecrafts, which obtained polarization maps for phase angles larger than a = 12°. The 
data show that the polarization in the B- and R-band for a ~ 90° reaches a level of about 
p ~ 50 % at the poles, while th e polarization is rather low (< 10 %) in the equatorial region 
(e.g. 



Smith and Tomaskdll984j ). 



For Saturn, iLyotI (Il929[ ) found a phase dependent polarization for the rings and some 
polarization for the atmosphere. Well established is the polarization of the ring system 
which shows for small phase angles < 7° a polarization of ab out p = —0.4 %, parallel to 



the s cattering plane, with some dependency on the phase angle (jjohnson et al.lll980l : iDoUfus 



19961 ). The disk of Saturn sho ws in the UV at 370 nm a strong radial limb polarization of 



more than 3 % near opposition (IHall and Rileylll974j ) as expected for Rayleigh scattering. In 



the visual the polarization is lower (^ 1 %), and predominately in N-S direction. The visual 
polarization shows significant temporal variations which may be seasonal. The poles show 
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DoUfus 
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Gisler et al. 


2003) 



Polari naetry of Saturn fo r large phase angl es ~ 30° — 150° was made w ith the spacecrafts 
Voyager 2 ( iWest et al.lll983l ) and Pioneer 11 (ITomasko and Dooselll984l ). These data show 
a strong wavelength dependence of the polarization with a low polarization (< 5 %) in the 
red, roughly 20 % in the blue and > 30 % (phase angle a = 68°) in the UV at 264 nm. A 
big step for ward is exp e cted f rom the imaging polarimetry of Saturn taken with the Cassini 
spacecraft. I West et al.l ( l2009l ) provide a first glimpse on a high quality Cassini polarization 
map of Saturn taken in 2003 for a phase angle of 61°. 

This paper presents imaging polarimetry and long slit spectropolarimetry of Jupiter 
and Saturn. In the next section a description of our observations and the data reduction are 
given. In Section 3 the polarization for Jupiter is described and analyzed while Saturn is 
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treated in Section 4. In Section 5 the observations are compared with model simulations and 
a qualitative interpretation is given for some polarization features. The results are discussed 
in the final section. 



2. Observations and data reduction 

Observations of Jupiter and Saturn were taken in November 2003 in spectropolari- 
metric mode with the ESO 3.6m telescope at La Silla, and in March 2002 and March 
2003 with polarimetric imaging at the McMath-Pierce solar telescope at Kitt Peak. Ob- 
servational parameters for J u piter and Saturn were taken or derived from data given in 



The Astronomical AlmanachI ( l2003l ) and they are summarized in Table [T] 



The illumination of the planet and the scattering geometry for the reflected light is 
defined by the phase angle a, which is the angle sun-planet-observer, and the position angle 
(PA) of the scattering plane 9 with respect to the central meridian (North-South direction) 
of the planet. The bright limb is on the East for 9 close to 90° and on the West for 9 close 
to 270°. 

The angles a and 9 are important parameters for polarimetry. The strength of polari- 
metric features depends on the phase angle a, and the orientation of the induced scattering 
polarization is often perpendicular or parallel to the orientation of the scattering plane 9. 
The scattering plane is essentially in East- West direction for 9 close to 90° or 270°. For the 
November 2003 observation of Saturn, the scattering plane is tilted by +12° with respect to 
the East- West direction. In this case a perpendicular or parallel polarization with respect 
to the scattering plane will produce besides a Q-polarization (in N-S orientation) also a sig- 
nificant [/-polarization component of f/ = 0.45(5. For a tilt angle of A^^ = —2° (e.g. Nov. 
2003 for Jupiter) this factor is f/ = -0.07(5. 

The apparent diameters ciN-s and ciE-w are used to convert locations x from the disk 
center (= sub-earth point) along the central meridian (CM) to radial distances tcm = x ■ 
c/n-s/S which can be converted to planetographic latitudes considering the ellipsoidal shape 
and the inclination of the planets (Tabled]). 



2.1. Spectropolarimetry 

Spectropolarimetric observations of Jupiter and Saturn were taken during the nights of 
November 29 and 30, 2003 with EF0SC2 at the ESO 3.6m telescope at La Silla. These data 
originate from the same run and instrument setup as the spectropolarimetry of Uranus and 
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Table 1: Parameters for the Jupiter and Saturn observations, a is the phase angle, 6 the 
orientation of the scattering plane, and tqu are distances from the sub-earth point on the 
central meridian. 



parameter 


Jupiter 




Saturn 




date (2003) 


March 9 


Nov. 29 


March 8 


Nov. 29 


observatory 


Kitt Peak 


La Silla 


Kitt Peak 


La Silla 


instrument 


ZIMPOL 


EF0SC2 


ZIMPOL 


EF0SC2 


a 


6.9° 


10.4° 


6.3° 


3.7° 


9 


269° 


88° 


272° 


102° 


polar axis inch 


+0.2° 


-1.5° 


+27.0° 


+25.0° 


lat. sub-earth point 


+0.2° 


-1.6° 


-32.0° 


-29.9° 


diameter (E-W) 


43.57" 


35.91" 


18.63" 


20.23" 


rcM limbs 


±20.38" 


±16.79" 


±8.59" 


±9.30" 


rcM south pole 


-20.37" 


-16.78" 


-7.49" 


-8.27" 


'"CM equator 


+0.1" 


-0.5" 


+4.2" 


+4.3" 


rcM ring inner edge N 






+6.45" 


+6.46" 


rcM ring outer edges 






±9.59" 


±9.61" 



Neptune from Joos and Schmid ( 2007al lbl). where descriptions of the measuring strategy and 



the data reduction are given. Here we provide only a brief outline and highlight some special 
points. 

EF0SC2 at the Cassegrain focus of the ESO 3.6m telescope is a multi-mode imager 
and grism spectrograph which can be equipped with a Wollaston prism and a rotatable 
super-achromatic half-wave plate for linear polarimetry and spectropolarimetry. A special 
slit mask can be placed in the focal plane which consists of a series of 19.7" long slitlets 
with a period of 42.2" appropriate to avoid overlapping of the spectra from the ordinary and 
extraordinary beam of the Wollaston prism. The width of the slitlets used for Jupiter and 
Saturn was 0.5". 

For extended sources the EF0SC2 instrument setup provides long slit spectropolarime- 
try and thus the intensity and polarization as function of wavelength A and position x along 
the slit. The orientation of the slit can be changed by rotating the whole instrument. The 
spatial resolution or the effective seeing of our data is about 1", as derived from the width 
of the spectra of the standard stars. 

Most of the Jupiter and Saturn data were taken with the slitlets along the central 
meridian. Some data were taken in East- West orientation (coordinates of the planet), but 
due to instrument flexure problems the calibration of the E-W slit data turned out to be less 
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accurate (see below). 

The grism (EF0SC2 grism^S) employed provided the wavelength range from 5200 to 
9350 A with a resolution of 6.4 A for a 0.5" wide slit. The data were recorded with a 2k x 2k 
CCD (ESO CCD#40) with a spatial scale of 0.157" and a spectral scale of 2.06 A per pixel. 
For A > 7000 A the CCD introduces an interference pattern. To smooth this pattern and 
to enhance the signal-to-noise ratio of the data the spectra are binned along the wavelength 
direction into 30 A bins. 



The linear polarization was measured in a standard way (e.g. iTinbergen and Rutten 



19921 ). with sets of four exposures taken with half-wave plate positions at 0°, 22.5°, 45° and 
67.5° respectively. For better data quality several sets were taken for each slit orientation. 
The exposure time per frame was 3 s for Jupiter, and 5 s for Saturn. Polarimetric standard 
stars were observed with the same instrumental setup to determine the polarization offset 
introduced by the telescope and to obtain the polarization angle zero point. The orienta- 
tion of the slit for standard star observations was always in celestial north-south direction. 
Exposures of a helium-argon lamp provided the wavelength calibration. 

The data reduction was performed with the MiDAS software package. For long slit 
spectropolarimetry it is important that the spectra of the ordinary and extraordinary beam 
are aligned with an accuracy of about 1/10 of a pixel in spatial direction, to ensure that no 
artificial polarization is introduced. This precision was achieved for the observations with 
the slit in N-S direction (N-S for t he planets was not far f rom N-S on the sky) using the 



standard star data as reference (see lJoos and Schmidll2007al ). The calibration of data taken 



with the slit in E-W direction suffered from instrument flectures and the standard stars 
spectra could not be used as alignment reference. Therefore the E-W spectra were aligned 
with respect to the edge of the slitlets. The achieved alignment precision is only about half 
as good as with the help of the standard stars. 

The instrumental polarization of EF0SC2 at the Cassegrain focus of the 3.6m telescope 
is low, and the polarimetric calibration of the data is straightforward. For the calibration we 
used the unpolarized standard star HD 14069 and the polarized standard star BD -|-25°727. 
The instrumental polarization was found to be less than 0.2 % in the central region of the 
field. The polarization angle calibration should be accurate to about ±2°. 

No flux calibration was attempted. Solar and telluric spectral features in the intensity 
spectra were corrected with the help of Mars observations taken with the same instrument 
configurati on. The ove r all slo pe of the intensity spectrum was adjusted to the albedo spec- 



trum from iKarkoschkal (119981 ) 



The observations of the linear polarization for Jupiter and Saturn are given as Stokes 
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parameters /(A, x), Q(A, x), U{\, x), where / is the total intensity and Q = Jq — /go, U = — 
1-4,5 are the polarized fluxes and Q/I, U/I the fractional (normalized) Stokes parameters. 
We also use the term radial polarization, e.g. Qr or Qr/I, indicating a polarization parallel 
(positive) or perpendicular (negative) to the slit through the disk center. The orientation of 
+Ur is rotated by 45° in counter-clock direction (North over East). 



2.2. ZIMPOL imaging polarimetry 

Imaging polarimetry of Jupiter and Saturn was taken in March 2002 and 2003 at the 
McMath-Pie rce solar telescope on Kitt Peak observatory, using ZIMPOL the Zurich Imaging 



Polarimeter (|Povellll995l : ICisler et all 120041 1. 



ZIMPOL is a high precision polarimeter, consisting of a fast polarization modulator, a 
polarizer, and a camera with a special masked CCD sensor. The polarization state of the 
incoming light is changed into a temporal polarization variation by the modulator, which is 
subsequently converted into an intensity variation by the linear polarizer. The masked CCD 
has periodically arranged open and covered rows. During integration the photo-charges are 
shifted back and forth between the open and one or more covered rows synchronously with 
the modulator. In this way the CCD camera has two or more image planes for the different 
polarization modes (e.g. Jq, /go ■■■), and all these images are registered with the same CCD 
pixels. Flat-fielding problems, differential aberrations, and alignment errors are minimized 
with this technique and the polarization can be determined with high sensitivity. 

This paper is focussed on the analysis of the spectropolarimetric signal of the strong 
polarization features of Jupiter and Saturn. The ZIMPOL imaging polarimetry is treated as 
useful auxiliary data because they cover the whole planetary disk of Jupiter and Saturn and 
Saturn's ring system in one exposure. This helps to verify and complete the spatial depen- 
dence of the spectropolarimetric data, for which one slit setting provides only incomplete 
spatial coverage. 

We have selected for our analysis the ZIMPOL imaging polarimetry of Jupiter and 
Saturn taken in March 2003 in a filter with a central wavelength of 7300 A and a band width 
of 200 A. Observational parameters for these data are given in Table [1] The pixel scale 
was 0.41 arcsec/pixel and the seeing was near 1.5 arcsec. In addition we use for Jupiter a 
ZIMPOL intensity image in the 6010 A "continuum" filter (width 180 A) for the calibration 
of the relative refiectivity scale on the disk. 

A drawback of the ZIMPOL / McMath-Pierce telescope observations are the position 
and wavelength dependent instrument polarization effects introduced by the inclined mirrors 
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of the solar telescope. We modelled the telescope polarization, but uncertainties at a level 
of Ag/J, Af/// ^ ±0.5 % remain. 

The instrument polarization can be corrected based on a priori assumption for the 
instrument and the targets. First we assumed that the instrumental polarization is field 
independent. Second we used the known polarization properties of Jupite r and Saturn. The 



reflected light shows r io circular polarization at a l evel \V/I\ < 0.1 % (jKemp et al.lll971 



Swedlund et al.lll972l : ISmith and WolstencroftI Il983l ). Further, the existing maps show an 
U / 1 polarization which is essentially anti-symmetric with respect to the central meridian an d 



a disk integrated [/-polarization close to zero \U / 1\ < 0.1 % (e.g. iHall and Rileylll968Ul974j ). 
This should also hold for the ZIMPOL data, because they were taken for an epoch where the 
scattering plane for Jupiter and Saturn was essentially in East- West direction; With these 
assumptions we can derive from our full Stokes polarimetry a cross talk corrected Q / /-image 
(see Fig. [1]). This Q//- image still includes all instrument polarization offsets induced by the 
inclined telescope mirrors. This offset was deri ved for our observations of Jupiter from the 
polarization phase curves of iMorozhenkd (119731 ). using the value Q/ 1 = —0.2 % for the disk 
center with an unce rtainty of less than ±0.1 %. For Sa turn the ring polarization is reasonably 



well known (see e.g. lDollfuslll996l : iJohnson et al.lll980l ) and for the offset correction we adopt 
a value of Q/I = —0.4 % for the polarization of the East and West cusps. 



3. Jupiter 

3.1. Imaging polarimetry 

Figure [1] shows the ZIMPOL Stokes-Q and U images of Jupiter taken in March 2003 in 
the 7300 A filter. A strong positive Q or U polarization flux is plotted white, while black 
indicates a strong negative polarization, and grey denotes little or no polarization flux. 

The Q flux image clearly shows the two strongly polarized poles with a polarization 
in N-S direction. Further there is a weak N-S polarization (parallel to the limb) near the 
equatorial limbs. In the center of the disk the polarization is slightly negative in East- West 
direction. The U image shows at the poles a positive or negative component on the East 
and West side of the central meridian. This indicates that the polarization at the poles is 
not in N-S but in radial direction. 

The fractional Q/I polarization at the poles reaches values of about ±7 % in the 7300 A 
filter. The disk-integrated (flux weighted) Q/Z-polarization is only ±0.2 %. The polariza- 
tion of the disk center is slightly negative according to the "i mposed polari z ation" offset 



calibration in accordance with the polarization phase curves of IMorozhenkd (119731 ). Our 



Fig. 1. — Q (left) and U (right) polarization flux images of Jupiter in a fllter centered at 
7300 A taken with ZIMPOL in March 2003 at Kitt Peak. North is up and East is left. 
The grey scale is normahzed to the central intensity and spans the range from —1.0 % 
(black) to +1.0 % (white). The lines in the Q-image indicate the sht positions for the 
spectropolarimetric observations with EF0SC2. 
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spectropolarimetric data (Sect. I3.3p confirm the calibration of tlie imaging polarimetry. 



3.2. Limb to limb profiles 



3.2.1. Comparison of North- South and East- West profiles 

Figure [2] shows the ZIMPOL 7300 A intensity and polarization profiles through the disk 
of Jupiter in North-South direction along the central meridian and in East- West direction 
along the equator. The plots give the radial polarization where positive Qr or Qr/I values 
stand for a polarization perpendicular to the limb. The position is indicated in arcsec with 
0" at the center of the apparent disk. The nominal limb position is at ±20.4" for the N-S 
profile and at ±21.8" for the E-W profile (see Tabled]). 

The N-S intensity profile in the 7300 A filter shows Jupiter's dark bands and bright zones 



struc t ure very si milar to many previous studies (e.g. IChanover et al.l Il996l : iMoreno et al. 
199ll : IWestlll979l ). The main feature of the equatorial intensity profile (dotted line) is the 
asymmetry due to the a = 6.9° illumination offset, with the bright limb in the West. 

The N-S and E-W fractional polarization profiles Qr/I illustrate the huge differences 
between the poles with a very high, positive limb polarization Qr/I ~ 8 % in the South 
and ~ 6 % in the North and the equatorial region with a small negative polarization. The 
polarization in the disk center is parallel to the scattering plane, which translates into a 
radial polarization Qr with a negative sign for the North-South profile and a positive sign 
for the East- West profile. 

The presented N-S polarization profile f or 730 A ag r ees well with ea rlier measurements 
for th e visual- red spectral region, e.g. from iLyotI (119291 ) . iDoUfusI (Il957l ). or iHall and Riley 
( 119681 ) . The polarization in the disk center depends on phase and it vari es from Q// 0.0 % 
for a = 0° to about Q/I ^ —0.5 % for a = 12° as described in detail in iMorozhenkd ( 119731 ). 



A negative limb polarization for the equator re gion was also previously reported for small 
phase angles < 6° and the red spectral region ( IDoUfusI Il957l : iGehrels et al.l 119691 ). In the 
UV/blue spectral region t he limb polarization at the equator is du e to Rayleigh scattering 
perpendicular to the limb (IHall and Rileylll968l : iGehrels et al.lll969l ). 



The radial polarization in f/r-direction is 5 to 10 times weaker and the signal obtained 
is dominated by systematic noise. The f/^-signal and the uncertainties are quantified for the 
spectropolarimetric observation. 



- 12 - 




'1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 'i-r-h-i--Ki'i I r 



o 5 

\ 

-i-> 
C 

o 



-25 -20 -15 -10 -5 5 10 15 20 25 

arcsec 

S (E) N (w) 

Fig. 2. — Observations of Jupiter from March 2003. Polarimetric profiles in N-S (full line) 
and E-W (dotted line) direction, taken in the 7300 A filter with ZIMPOL. The top panel 
shows the intensity /, the middle panel the fractional (radial) polarization Qr/I, and the 
bottom panel the polarization flux Qr- 
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3.2.2. Polarization measurements for the central meridian 

The long-slit spectropolarimetric measurements provide polarimetrically well calibrated 
profiles for the central meridian for all wavelengths between 5300 and 9300 A. 

Figure [3] shows profiles for the continuum at 6000 A, spectrally averaged from 5900 to 
6100 A, (solid line) and the deep methane absorption band at 8870 A spectrally averaged 
from 8800 to 8940 A. Spectral averaging was done to enhance the signal to noise. 

The intensity profiles show belt and zones at 6000 A and a significant li mb brightening 



in the 8870 A methane ba nd in good agreement with previous studies (e.g. Moreno et al. 



199ll : IChanover et al.lll996r ). 



The Qr/I and Qr polarization profiles from November 2003 in Fig. |3]are very similar 
to the profile from March 2003 shown in Fig. [2] with Qr/I ~ 10 % at the south pole for 
both spectral bands. At the north pole Qr/I for 6000 A reaches a maximum value just 
above 8 %, while the maximum for 8870 A is slightly less than 6 %. It is interesting to note 
that the South pole shows not only a stronger polarization than the North pole, but also a 
stronger limb brightening in the methane absorption A8870. The polarisation is negative in 
the center of the disk. The sign change occurs at about ±12.5 arcsec, corresponding to a 
Jovian latitude of about ±59°. 

It is not trivial to quantify accurately the measured polarization signal. The peak 
polarization at the limb depends significantly on the spatial resolution (or the effective seeing) 
and it is difficult to compare our measurements with previous data which often depend on 
the spatial resolution of the measurement. 

In order to provide quantitative results we split the profile into the polar sections with 
positive polarization S± and N± and sections with negative polarization S— and N— . The 
section borders are set at ±12.5 arcsec, where the signs change, and at ±5 arcsec, the approx- 
imate end of the slitlets. For measurements of the central region without spectropolarimetric 
coverage we interpolated the curves base d on the ZIMPOL A7 300 imaging polarimetry and 



intensity profiles from the literature (e.g. IChanover et al.lll996l ). 



Table [2] gives the fiux /, polarization fiuxes Qr and Ur, and the fiux weighted fractional 
polarization Qr/I and Ur/I for the wavelengths 6000 A, 8200 A, and 8870 A. The fiux is 
given as ratio I/Isia relative to the fiux for the entire slit. This quantity does not depend on 
uncertainties in the absolute fiux calibration. 

The values Qr/I and Ur/I for the S± and N± regions in Table [2] depend strongly on 
the exact size of the integration interval due to the steep intensity slope I{x) at x = ±12.5". 
Much less critical is the determination of the total polarization fiux Qr from x = ±12.5" 
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Fig. 3. — Observations of Jupiter from November 2003. Intensity J, fractional polarization 
Qr/1 1 and polarization flux Qr profiles through the central meridian (N-S) of Jupiter for the 
continuum at 6000 A (full line) and the deep methane absorption band at 8870 A (dotted 
line). The / and flux for the methane band are multiplied by 10 for visibility reasons. 
At the top the five selected spatial regions are indicated. 
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Table 2: Fractional intensity ///sut, radial polarization Qr/I and Ur/I for Jupiter for the 
sections S+, S— , center, N— , N+ along the central meridian with boundaries as defined in 
the first tw o rows . Ag is the geometric albedo for the considered wavelength range from 



Karkoschkal (119981 ) and / the derived refiectivity (see text). 





s+ 


S- 


center 


N- 


N+ 


total slit 


location 














•^min 


-20" 


-12.5" 


-5" 


+5" 


+12.5" 


-20" 




-12.5" 


-5" 


+5" 


+ 12.5" 


+20" 


+20" 


Ax/xslit 


0.128* 


0.223 


0.298 


0.223 


0.128* 


1.000* 


continuum 5900-6100 i 


^, A = 


0.54 








///slit 


0.067 


0.250 


0.366 


0.250 


0.067 


1.000 


f 


0.31 


0.66 


0.72 


0.66 


0.31 


0.59 


Qr/I [%] 


+4.36 


-0.56 


-0.44 


-0.44 


+3.95 


+0.14 


Ur/I [%] 


+0.35 


0.00 


-0.01 


-0.05 


-0.31 


-0.02 


Qr/Ism[%] 


+0.292 








+0.265 




continuum 8100-8300 A, A„ = 


0.48 








///slit 


0.057 


0.246 


0.368 


0.255 


0.074 


1.000 


f 


0.23 


0.57 


0.64 


0.59 


0.30 


0.52 


Qr/I [%] 


+3.31 


-0.68 


-0.42 


-0.61 


+1.73 


-0.16 


Ur/I [%] 


-0.40 


-0.02 


+0.01 


+0.02 


-0.24 


+0.01 


g.//siit[%] 


+0.189 








+0.128 




CH4-band 8 


800-8940 A, A„ = 


0.05 








///slit 


0.095 


0.203 


0.388 


0.229 


0.096 


1.000 


f 


0.039 


0.047 


0.068 


0.053 


0.039 


0.052 


Qr/I [%] 


+5.69 


-1.16 


-1.07 


-1.21 


+2.79 


-0.15 


Ur/I [%] 


+0.61 


+0.07 


+0.26 


+0.02 


-0.63 


+0.10 


Qr//slit[%] 


+0.541 








+0.268 





*: Ax/xsiit considers only the slit section located within the limb at ±16.79". 
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(where Qr ~ 0) to the hmb. The measured value can be related to /sut, giving Qr{S+) j Iswt 
and Qt{,N+) / Iswt according to Qr/Isiit = I/hnt ■ Qr/I (also for Ur/Isia)- This relative limb 
polarization flux is an interesting quantity for long term variability studies, since it does 
not depend much on the spatial resolution of the observations and can be derived from the 
counts data without flux calibration or conversion into normalized reflectivity. It remains to 
be determined whether Qr(5'+)//siit and Qr{N+)/Is\it depend on phase angle. 

A conversion of the fractional intensity I/Isut into reflectivity / can be made. If the 
average reflectivity along the central meridian (/sut) is known, then the reflectivity in a bin 
is 

e _ / f \ -^/-^slit 

J — \Js\\t) — -, • 

X/ XsWt 

From the full disk image we can derive the ratio A = (/disk) / (/sUt) between the average 
reflectivity for the full planetary disk ( = the geometric albedo Ag) and fswt- With the 
geometric albedo from the literature one can determine 

(/slit) = • 

We derive A = 0.92 from the ZIMPOL image for the "continuum" filter centered at 6010 A 
(width 180 A). This value is not m uch different from ALam = 0.85 for a perfectly white 



Lambert sphere. iKarkoschkal (Il998l ) gives for this wavelength a geometric albedo of Ag = 
0.59 for data taken in 1995. This value can be used for our calibration because the global 
reflectivity variations of Jupiter are small (< 5 %) and also the phase dependence of the 
reflectivity (~ 1.5 % for a 0° — 10°) can be neglected. A = 0.92 can also be employed to 
derive the reflectivities for the continuum wavelength 8200 A, because our intensity profiles 
along the central meridian have a very similar shape for wavelengths from 5300 A to 8700 A. 

More difficult is the calibration of the reflectivity for the strong methane band at A = 
8870 A. This profile shows a relatively high reflectivity at the equator and the poles and 
relatively dark mid- latitudes (see Fig. [3]). The overall reflectivity profile is fiatter than for 
the continuum wavelengths 6010 A or 8200 A, but not completely fiat (A = 1) as for a 
refiecting disk. Strongly absorbing atmosphere s, like Jupiter in a strong CH4 band, have 



a rather constant refiectivity over the disk (see iBuenzli and SchmidI 120091 ). Based on these 



considerations we adopt A = 0.96 for the 8870 A-methane band. 

The uncertainties for the fractional polarization values Qr/I and Ur/I given in Table 
|2] are mainly due to systematic errors like instrument calibration or inaccuracies in the slit 
positioning. The derived Qr/I and Ur/I values should be accurate to A{Qr/I), A{Ur/I) ~ 
±0.1 % to 0.2 %, except for the Ur/I-valnes for S+ and N+, for which a Q ^ U cross-talk 
error or/and slit positioning error at a level of up to Af/,.// ^ ±0.3 % to ±0.6 % could 
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be present. We therefore conclude that our data show no significant Ur component for the 
central meridian of Jupiter. 



3.3. Spectropolarimetry 



The EF0SC2 observations provide the spectropolarimetric signal for each point along 
the slit. A general picture is provided in Fig. H] which shows averages for the spatial regions 
S+, S— , N— , and N+, defined in the previous section. 

Th e refiection spectra /(A) are color-calibrated with respect to the full disk albedo spec- 
trum of iKarkoschka (119981 ). Spectroscopic features are very similar for all regions (see also 
Cochran et al.lll98l[ ) but systematic differences do exist. For example the equivalent width 
of the strong A8870 CH4 band is smaller for the limb regions 8+ and N+, when compared 
to mid latitudes S— or N— . This is just another manifestation of the limb brightening effect 
for CH4 A8870 seen in Fig. El 

The fractional polarization Qr/I{X) shows for the South polar region S+ a decrease of 
the continuum polarization with wavelength from about Qr/I = 4.6 % at 5300 A to 2.8 % 
at 9300 A. The decrease is much steeper for the northern limb from Qr/I = 5.0 % to 1.2 %. 
In the strong absorption bands Qr/I is enhanced with respect to the adjacent continuum, 
most prominently in the CH4-band at A8870. 

The fractional polarization for the mid-latitude regions S— and N— is low, but negative, 
and the continuum polarization increases slightly from about Qr/I = —0.3 % at 5300 A to 
—0.6 % at 9300 A. The polarization in the methane bands is also enhanced, in this case more 
negative (Fig. SI Table [2]). 

The polarization flux spectrum Qr(A) = I-Qr/I shows a reduced signal at the position of 
the strong absorption bands. This indicates that the enhancement of Qr/I in the absorption 
bands is smaller than the reduction in /, so that the polarization flux spectra still exhibit a 
reduction in |(5(A)| at the wavelengths of the absorptions. 

Essentially no previous spectropolarimetric data fr om Jupiter is a vailab le in the litera- 
ture apart from the multi-fllter aperture polarimetry of iGehrels et al. fll969h and the detec- 
tionof a weak differential polarimetric signal due to a methane ban d by lSmith and Wolstencroft 
(Il983l ) measured with full disk observations. iGehrels et al.l (Il969[ ) found a similar wavelength 
dependence for the polar limb polarization. Speciflcally, in 1960/63 they measured for longer 
wavelength (> 5000 A) a higher polarization at the northern than at the southern limb - 
the opposite N-S asymmetry when compared to our observations. They also found that the 
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polarization asymmetry is reversed for shorter wavelength (< 5000 A), and the South pole 
shows a higher polarization. Such a reversal of the N-S polarization asymmetry also seems 
to be present at the short wavelength end of our data. 

3.3.1. Spectropolarimetry for the "extreme" polar limbs 

The S+ and N+ spectropolarimetry provides only a very coarse characterization of the 
"average" limb polarization. For this reason we explore in this section the spectropolarimetric 
signal at the "extreme" polar limbs, where we find the maximum of the fractional polarization 
Qr/I. We select a small spatial bin from ±16" to ±17", corresponding to planetographic 
latitudes from about ±73° out to the polar limbs (see Fig. [5]). 

For the South pole the fractional polarization reaches Qr/I = 9.6 % for the continuum 
m the 5300 to 6100 A region and even 10.0 % in the A6190 CH4-band. Previous studies 
report maximum fractional polarizations of not more than Qr/I = 7 — 8 %, most likely 
because the spatial resolution was worse than the 1" achieved with our observations. This 
illustrates the problem of the spatial averaging. The overall spectral slope of the fractional 
polarization spectra Qr/I{X) at the extreme polar limb is flatter for the South than for the 
North, similar to the S± and N± regions. Interestingly the polarization enhancement in the 
strong methane bands is much less pronounced in the extreme polar limb data. 

3.3.2. Spectropolarimetric signal at the equatorial limb 

The spectropolarimetric observations of Jupiter taken with a slit in East- West direction 
cover the regions from the eastern limb at —18.0" to —12.8" and from +10.7" to the western 
limb at ±18.0" (Fig. [I]). 

For both limbs the fractional polarization spectra show a positive polarization Qr/I ~ 
0.5 % for the shortest wavelengths, and Qr/I ^ 0.0 % ±0.2% for A > 7000 A, except for the 
strong CH4 band at 8870 A, where the polarization is approximately Qr/I ~ ±0.4 %. 

The spatial profiles show a general radial decrease in the fractional polarization from 
slightly positive values, Qr/I ~ ±0.5 % to 0.0 %, at the inner edges of the slitlets (at —12.8" 
and ±10.7") towards zero or slightly negative values, Qr/I ~ 0.0 % to —0.5 %, further out at 
the eastern limb and out to ±16" on the western side. For a given distance from the center, 
the polarization is slightly more negative {Ap ^ —0.2 %) on the eastern side. 

This behavior is qualitatively in agreement with the E-W polarization profile derived 



0.0 




J I I I I I I I I I I I I I I I I I I L 

6000 7000 8000 9000 

wavelength [A] 



Fig. 4. — Spectropolarimetry of Jupiter for the S+, S— (solid lines) and the N— and N+ 
(dashed hnes) regions as defined in Fig. [3l The intensity /(A) for the polar region S+ and 
N+ are multiplied by a factor of 3 with respect to S— and N— for visibility reasons. The 
middle panel gives Qr/I{\) and the bottom panel the polarization flux Qr- 
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Fig. 5. — Fractional polarization spectra Qr/I and Ur/I for the "extreme" southern (solid) 
and the northern (dashed) limb. 
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for the ZIMPOL polarimetry for March 2003. The E-W asymmetry is switched most hkely 
because the bright hmb and terminator have switched sides between March 2003 and Nov. 
2003. Also the small change in the overall polarization level for the phase angle a = 6.9° in 
March 2003 (less E-W polarization in the equator re gion) and a = 10.4 ° in Nov 2003 (more 
E-W polarization) is similar to previous studies (e.g. iMorozhenkd Il973l ) . 



4. Saturn 

4.1. Imaging polarimetry 

The appearance of Saturn depends strongly on the inclination of the planet and its ring 
system. In 2003 the inclination of the polar axis with respect to the celestial plane was 
27°, close to the maximum inclination possible for Earth-bound observations. An EF0SC2 
acquisition image shows Saturn for the November 2003 run (Fig. [6]). For this inclination the 
South pole is well visible about 1" inward from the limb (see Tabled]). The equator is almost 
half way between disk center and northern limb while the ring system covers the latitudes 
northward of +20°. 

Figure [7] shows a Stokes Q-polarization flux image in the 7300 A filter taken with 
ZIMPOL in March 2003. One can recognize the negative (E-W) polarization of the ring, 
a weak positive feature at the southern pole, a polarization flux close to zero in the south, 
corresponding to latitudes around —60°, and a dominant strong positive feature at the 
equator. 

The [/-image corresponding to Fig. [7] shows near the southern limb negative ?7// ~ 
—0.3 % and positive ~ +0.3 % features on the east and west sides of the South pole, 
respectively very much like the [/-pattern for the poles of Jupiter. This indicates that the 
polarization features at the South pole have a radial polarization direction. For the rings of 
Saturn there is essentially no [/-signal visible. 

The equatorial polarization reaches in the 7300 A filter a maximum fractional polariza- 
tion of Q/I ~ +1.6 % and extends about 4" to 5" in east-west direction. This equatorial 



featu r e is remarkab le since it was not present in early measurements (e.g. iHall and Riley 



1974 : IPoUfud Il996 ). A similar signal is also present in the 6010 A imaging polarimetry 



from the sa me date and in the 5500 A, 6010 A, and 7300 A filter ZIMPOL polarimetry from 



March 2002 (iGisler et al.ll2003l ). Also our spectropolarimetric observations of November 2003 



confirms the presence of this transient signal. 



The [/-signal at the equator is weak with a fractional polarization of [/// ~ —0.3 % 
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at the east limb, —0.1 % at the meridian, and about +0.1 % on the west hmb. These are 
only small signals when compared to the fractional ^//-polarization of ~ +1.5 % at the 
intersection of meridian and equator. This indicates that the equatorial polarization extends 
from about —30° to +30° in longitude along the equatorial band and has a predominant 
orientation in north-south direction. 



4.2. Profiles for the planet and the ring 

4-2.1. Comparison of North- South and East-West profiles 

Figure [8] shows the North-South and East- West profiles extracted for the 7300 A CH4 
filter ZIMPOL observations of March 2003. The middle panel shows the fractional radial 
polarization Q,/!- The N-S profile for Qr/I is positive everywhere on the planetary disk 
and negative for the ring. There are two peaks, one of Qr/I ~ +1.0 % at the southern 
pole (at approximately —8.5") and a stronger peak of Qrll ~ +1.6 % at the equator (+4") 
coincident with the intensity maximum. For the E-W direction the radial polarization is 
positive for the ring (parallel to the slit) and negative for the planet. The strongly negative 
Qrll polarization between planet and ring is probably a spurious effect because in these 
gaps the photon statistics is low and the impact of scattered light, and residual systematic 
noise effects are large. The polarization flux profile Qr shows that the equatorial polarization 
feature is really strong. 

As described in Sect. 12.21 the calibration for the ZIMPOL-polarimetry is based on the 
assumption that the fractional ^//-polarization is —0.4 % for the ring. The good agreement 
between the ZIMPOL and EF0SC2 data support this calibration procedure. 

4-2.2. Polarization measurements for the central meridian 

Figure [9] compares the meridional profiles for the strong methane band A8870 and the 
continuum at A6000. For the A8870 CH4 band the fractional polarization Qr/I is positive 
for the southern latitudes and negative at the equator and the ring. It appears that Qr/I 
has a maximum at the southern pole at about —7.5". The polarization flux Qr for A8870 
illustrates the weakness of the polarization signal at this wavelength. In Qr no polarization 
flux maximum is visible at the south pole. Interestingly there is a small negative polarization 
for the CII4 wavelength near the equator in contrast to the strong positive polarization signal 
at shorter wavelengths. 



Fig. 6. — EF0SC2 acquisition image of Saturn taken on Nov. 30, 2003. North is up and 
East is to the left. The hmb of the planetary disk, as well as the position of the south pole 
and the limits of the A and B rings are indicated. In addition, the two slit positions for 
spectropolarimetry are marked with white rectangles which overlap at the center. 




Fig. 7. — ZIMPOL Stokes Q-flux image of Saturn, obtained in the methane 730 nm filter. 
In the south the planet occults the ring, except for a small rim with a width of about 1.5". 
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Fig. 8. — Saturn N-S (full line) and E-W profiles (dotted line) for the intensity (top), the 
fractional radial polarization Qr/I ■, and the Qr polarization flux for the wavelength A7300A. 
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Fig. 9. — Polarimetric profiles of Saturn along the central meridian for the continuum wave- 
length region 6000 A (full line) and the deep methane absorption band at 8870 A (dotted 
line) . 
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We select four regions along the the central meridian with a width of 2" as indicated 
in Fig. [9l These regions are "S" for southern limb, "C" for the center of the disk, "E" for 
the equator, and "R" for the ring and their relative intensity and fractional polarization are 
given in Table El 

For the interpretation of the polarization in Table |3] the orientation of the scattering 
plane has to be taken into account. For the November 2003 observations the scattering 
plane is tilted by +12° with respect to the E-W direction of Saturn. Therefore a scattering 
polarization perpendicular to the scattering plane produces besides a Q-signal also an U 
signal at the level of U = 0.45 Q. For example, the A6000 polarization for the "S" region 
has an orientation of dp = +5°, compatible with 6p = 0° for a polarization in N-S direction 
(or radial), while the strong polarization in the equatorial region "E" has an orientation of 
6p = +11°, perpendicular to = 102° of the scattering plane. 

For the polarization of the E and W cusps of the ring the litera ture gives values of 



0.4 % parallel t o the scattering plane for phase angles near 4° (e.g. lDollfuslll979l . Il996 



Johnson et al.lll980[ ). For the ring "R" we measure for A6000 and A8100 a. Q/I polarization 
component between —0.1 % and —0.2 %, or a polarization which is low er (less negative) by 
about 0.2 % than expected. This effect has been previously described by lKemp and Murphy 



(Il973l ). and at that time it was explained by an admixture of +(5//-polarized light from 
Saturn transmitted through the ring. Since we know now that the optical depth of the ring 
is rather high, we attribute this decrease in polarization to light from the adjacent equatorial 
region scattered by the instrument. In the dark CH4 band the reflection from the planet 
is strongly reduced, so that the ring polarization is less diluted by scattered light from the 
planet. 

The uncertainties for the fractional polarization due to instrument calibration and slit 
positioning errors are about AQr/I,AUr/I ~ ±0.2 %. Photon noise errors are negligible 
compared to calibration errors and unidentified instrumental effects. 

Our relative I/Isnt for Saturn can be co nverted into reflectivities /a using the HST ob- 
servations from December 2002, published by iPerez-Hoyos et al.l (120051 ). From the electronic 
figure we obtain f^j^o = 0.72 + 0.01 for the sub-Earth point "C". Temporal changes in fx for 
Saturn can be significant, due to the strong inclination of the planet and the variable shadow- 
ing by the ring system. But the data from December 2002 are useful for us because no strong 
reflectivity changes were noticed between D ec. 2002 and Aug. 2003 (IPerez-Hovos et al.ll2005r) 
and between March 2003 and March 2004 (IKarkoschka and Tomaskdl2005l ). 



Reflectivities for other wavelengths are deduced with the global albedo spectrum of 
Karkoschkal fll998h . according to fx = /6750 ^s(A)/v4g(6750), using Ag{6750) = 0.57. The 
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Table 3: Relative intensity ///gUt, reflectivity /, and fractional radial polarization Qr/ 1 and 
Ur/I for Saturn for the four 2" wide regions S, C, E, R on the central meridian and the 
entire meridian. The first row gives the center of these meridional sections and the second 
row the corresponding planetographic latitude. 







s 


C 


E 


R total slit 














location 




•2'cent 




-9.3" 


0.0" 


+4.0" 


+8.0" 






latitude 


-65°(") 


-30° 


0° 


+30° 












continuum 


5900-6100 


A, Ag — 


0.54 


///slit 




0.037 


0.124 


0.129 


0.084 


1.000 




/ 






0.68 


0.71 


0.46 






Qr/I [ 


%] 


+ 1.40 


+0.34 


+ 1.65 


-0.10 


+0.54 




Ur/I [' 


%] 


+0.24 


+0.08 


+0.68 


-0.17 


+0.16 










continuum 


8100-8300 


A, Ag — 


0.56 


///slit 




0.027 


0.125 


0.130 


0.088 


1.000 




/ 






0.71 


0.74 


0.50 






Qr/I [ 


%] 


+0.40 


+0.16 


+0.02 


-0.21 


+0.05 




Ur/I [' 


%] 


+0.15 


+0.02 


-0.22 


-0.16 


-0.05 












CH4 88( 


DO-8940 A, 


Ag = 0.07 


///slit 




0.045 


0.047 


0.102 


0.301 


1.000 




/ 






0.09 


0.20 


0.6 






Qr/I[ 


%] 


+0.27 


+0.19 


-0.35 


-0.40 


-0.18 




Ur/I [ 


%] 


+0.36 


-0.04 


-0.74 


-0.50 


-0.38 





a: southern limb beyond the south pole at longitude +180° from the central meridian 
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reflectivities for tlie slit sections S, E, and R scale like the I/Isia values (Table |3]). 



4-2.3. Spectropolarimetry for Saturn 

The spectropolarimetric signal for the southern limb (S), disk center (C), equator (E), 
and the ring (R) are plotted in Fig. [TUJ The spectra are averages for 2" wide spatial regions 
as indicated in Fig. [91 

The intensity spectra S, C, and E from the planetary disk show all the well known 
CH4-bands. The ring spectrum R is essentially featureless apart from very weak features at 
the wavelength of strong CH4 bands, which can be attributed to scattered light from the 
adjacent planetary disk. 

A strong fractional polarization {Qr/I > 1 %) is present at the south pole, "S", and 
equator, "E", respectively. Qr/I decreases steeply towards longer wavelengths from about 
Qr/I ^ +2.2 % at 5300 A to about +0.2 % at 8000 A, and is close to zero for longer 
wavelengths. The fast drop in polarization in Qr/I or in the polarization flux Qr with 
wavelength i s typical for a Rayleigh sc attering layer above a weakly or non-polarizing cloud 



surface (e.g. iBuenzli and Schmidll2009l ) 



Qr/I for the disk center is low, only +0.4 % at 5300 A, and it decreases steadily towards 
zero. It is not clear, whether the negative feature i n the CH4-band at A8870 is real or just 



a spurious effect. As shown in previous studies (e.g. iKemp and Murphy! 1 19731 : lDollfuslll996l ) 
the fractional polarization of the ring depends only little on wavelength. A weak feature is 
visible in the strong CH4-band A8870, which could be due to enhanced contamination from 
scattered light. 

For the "E" and "S" regions there are essentially no narrow spectral features visible 
in Qr/I- In particular there is no clear enhancement in strong methane bands like A 72 70 



or A88 70, as observed for Jupiter (see Sect. 13. 3p . or Uranus and Neptune (jjoos and Schmid 



2007af l. 



5. Investigation of the strongest polarization features 

Our polarimetry of Jupiter and Saturn has revealed two surprising polarization features: 
a very high limb polarization reaching a maximum of more than 9 % between 5300 A and 
6500 A for the South pole of Jupiter, and a strong equatorial polarization for Saturn. We 
explore whether these observational results are compatible with simple scattering models. 
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Fig. 10. — Spectropolarimetry for the four regions of Saturn: south pole (solid line), center 
(dashed line), equator (dot-dashed line), and ring (dotted line). Top panel: intensity spectra 
/(A) (counts). Middle panel: fractional polarization Qr/I{X). Positive is parallel to the slit 
and negative is perpendicular. Bottom panel: polarized flux Qr{)^)- 
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For a detailed characterization of the scattering particles in Jupiter and Saturn extensive 
modelling would be required. Unfortunately there exist up to now only few limb polarization 
models - essentially only the grid for Rayleigh scattering atmospheres bvlBuenzli and Schmid 
(I2OO9I ) and a few previous, mostly analytic results as summarized in ISchmid et al.l (l2006al ). 
Therefore it is not well explored yet how the limb polarization depends on the scattering 
phase matrix of the haze particles, the stratification of the atmosphere, and the optical 
depth of absorbers. Therefore, our model fitting remains ambiguous without extensive model 
simulations which are beyond the scope of this paper. 



5.1. Polarization model for the poles of Jupiter 



Model calculations for Jupiter were carried out in order to explain the observed peak 
limb polarization of more than 9 % in R-band and more than 9.5 % in the V-band at 
the South pole. Considering that the seeing degrades this polarization, the maximum limb 
polarization must be well above 10 % . Rayleigh scattering models yield up to 10 % limb 
polarization (e.g. ISchmid et al.l |2006aJ) but only for highly absorbing atmosphere models. 



which are not appropriate for the reflected intensity seen on the poles of Jupiter. 



Detail ed scattering r no dels for the polarization of Jupiter were presented by lSmith and Tomasko 



( Il984j ) and iBraak et al.l ( 120021 ) . but only for r nid and low latitudes . No detailed scattering 
models exist for the polarization at the poles. ISmith and Tomaskol ( 1l984j ) made a simple fit 
to the polarization measured in the red with the Pioneer spacecraft for a Rayleigh scattering 
layer with single-scattering albedo of w = 0.983, optical thickness r = 0.5, and a surface 
albedo of As = 0.67. However, this Rayleigh scattering model yields a maximum limb po- 
larization of only 7.3 %, or fa 6.5 % if the degradation by the seeing is considered, whereas 
our measurements show a much higher fractional polarization. 

We calculate the polarization along the central meridian of Jupiter, considering three 
zones: the polar S+ and N+ zones, and a central zone, corresponding roughly to the S— , 
N— , and center regions in Table |21 



m 



For the S-f and N -| - zon es we use the Monte Carlo multiple scattering code described 



Buenzli and SchmidI (120091). The chosen atm osphere structure is very similar to the haze 



model presented by lSmith and Tomaskol (119841 ) for low latitudes. They determined haze and 
gas properties for the South Tropical Zone and South Equatorial Band with an atmosphere 
consisting of a top gas layer Gl, a scattering haze layer H, a lower gas layer G2, and an 
optically thick surface layer S at the bottom. We fit the polar regions S+ and N+ with this 
model. For the wavelength 6000 A the optical depths for the two gas layers are tqi = 0.011 



- 31 - 



and tg2 = 0.018, and the single scattering albedo is lug = 0.976, calculated for a methane 
abundance of 0.18 %. Thus, the haze layer is geometrically thin and located at a pressure 
level of 290 mbar, while the opaque surface layer is at 760 mbar. The continuum polarization 
at 6000 A does not depend much on the exact height of the haze and cloud surface layers 
and e.g. 40 mbar or 1200 mbar, respectively, would not notably change the intensity and 
polarization. The situation is different for the polarization in the CH4 bands which depends 
in various ways on the pressure level of the different layers. 

The free parameters of our model are: the albedo cloud layer As assumed to be a grey 
Lambert surface, the optical thickness of the haze layer Th, and the haze parameters, which 
are single scattering albedo cuh, single term Henyey-Greens tein asymmetry pa rameter gi^, 



and maximum polarization for right angle scattering (see iBraak et al.ll2002l ). The polar 
model is independent of latitude, but the incidence and viewing angles produce a latitude 
dependence in the reflected polarization and intensity. 

For the reflected intensity at lower latitudes the same model is used, but for the frac- 
tional polarization an "artificial" constant value of Qr/I = —0.7 % is adopted. The negative 
polarization is introduced to fit the transition between the positively polarized polar zones 
and the negative central zone. The borders and are free parameters which are deter- 
mined in the data fitting process. 

In order to describe the smearing of the signal due to atmospheric seeing a nd instrumen- 



tal light scattering, the "discrete" three zone mode l is convolved with a Moffat (jMoffatlll969[ ) 



point spread function (PSF). In the formalism of iTrujillo et al.l (120011 ) this PSF includes a 
/9-parameter which describes the scattering wing. A small /3 implies strong wings, a Gaus- 
sian is obtained for /3 — )■ 00, while atmospheric turbulence theory predicts /3 = 4.76. For our 
observations we derive (3 = 1 from the residual light outside the nominal limb, indicating 
significant scattering in the instrument. 

In Fig. [11] the observed intensity and polarization profiles for the central meridian for 
6000 A are compared with the model fit. At the poles |r| > 0.8 the match is satisfactory 
except for the fractional polarization outside the limb \r\ > 1, where the statistical errors are 
large because / — >^ 0. At low latitudes there are some discrepancies because we did not try 
to fit the band structure for the intensity, and we adopted a constant value for the fractional 
polarization. 

A good fit for the limb polarization at A = 6000 A for both poles is obtained for strongly 
polarizing ph = 1 haze with a low absorption = 0.99 and asymmetry parameter g-^ = 0.6. 
Such scattering paran ieters are typical for ag gregate particles as proposed to be present in 



Jupiter and Titan by IWest and Smith! (jl99ll ). Small Mie scatterers with diameters much 
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smaller than the wavelength of the scattered light have similar scattering parameters. But 
the scattering cross section for small spheres is much higher in the blue and one would expect 
a fractional limb polarization which decreases rapidly with wavelength (as seen for the pole 
and the equator of Saturn). The scattering cross section and therefore the scattering optical 
depths of the haze layer are expected to decrease only slowly with wavelength for aggregate 
particles and they can therefore explain the rather gentle decrease in the Q/ 1 and Q-spectra 
towards longer wavelengths observed for the poles of Jupiter (see Fig. H]). 

The measured limb polarization at the poles for r = ±0.96 is about 9.8 % in the South 
and 8.4 % in the North. The North-South differences can be explained by different optical 
depths T\i{N+) = 0.72, rh(5'+) = 1.1 for the haze layers. 

Our model fit includes the smearing due to seeing and the modelling indicates that 
the intrinsic limb polarization reaches maxima of about 11.5 % in the South and 10.0 % 
in the North. Such a high limb polarization is probably only possible with particles having 
a scattering phase function with reduced backscattering, with an asymmetry parameter 
comparable to ~ 0.6 as used in our model. 

Our model fit is not unique. However, many parameters are already well constrained. 
From the parameter space explored by us it seems very likely that the asymmetry parameter 
lies in the range 0.8 < (^h < 0.5 and the single scattering polarization is ph > 0.8. 

The locations of the transitions and between the highly polarized poles and the 
negatively polarized central zone are at r = ±0.82 ± 0.01, identical for both hemispheres. 
This correspond to a planetographic latitude of ±58° ±3°. In the South the transition 
between the S± and the central region is compatible with an unresolved discontinuity. In 
the North the transition is more gradual, with a weak wing of positive polarization toward s 



lower latitude. This is in qualitative agreement with Fig. 14 in ISmith and Tomaskd (119841 ). 
who measured the latitudinal polarization dependence with Pioneer for large phase angles 
(82° and 98°). Their red filter data show a steep increase in the polarization in the South 
from about 10 % at a latitude of —55° to about 33 % at —65° and a more gradual increase 
from 13 % at +55° to about 25 % at +65° in the North. 

The surface albedo As = 0.75 is also quite well constrained, since a high value As > 0.9 
would significantly reduce the resulting fractional polarization and a low value As < 0.6 
would underpredict the reflectivity at the poles. It also seems quite safe to explain the 
North-South asymmetry in the polar polarization with a difference in the optical thickness 
of the polarizing haze layer. 



Finally, it is interesting to n ote that the haze model for low latitudes with ooh = 0.95, 
Ph = 0.9 and Qh = 0.75 from ISmith and Tomaskd (119841 ). which was also adopted by 
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Fig. 11. — Comparison of the model fit (dashed hne) with the observed 6000 A N-S profile 
(dotted line) of Jupiter. 
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Braak et al.l (120021 ) . cannot fit the polarization at the poles. This points to distinct dif- 
ferences between the haze particles at the poles and at lower latitudes. 



5.2. The transient polarization feature at the equator of Saturn 

For Saturn we investigate the nature of the surprisingly strong transient polarization 
feature at the equator Q^//(6000 A) = +1.65 % (Table [3]), present in our 2003 (and 2002) 
data. The transient nature of the feature could be related to the high inclination of ~ 26° of 
the Saturn system during our observations. In the following we explore possible explanations 
for this feature. 

At exact opposition no polarization signal is expected from the center of the planetary 
disk for symmetry reasons. Our spectropolarimetric observations were taken at phase angle 
a = 3.7°. The equatorial polarization cannot be ascribed to a phase effect because for 
this phase one expects only a negligible positive polarization for small, positive polarizing 
scattering particles, or a small Q/I ~ —0.2 % negative polarization for e.g. large negatively 
polarizing scattering particles. 

Higher order scatterings produce a radial limb polarization, but this initially increases 
only slowly when going from the center of the disk towards the limb. A semi-infinite, con- 
servative {u = 1) Rayleigh scattering atmosphere p roduces a fractional limb polarization of 
only O r/ 1 = +1.3 % at r/rumh ~ 0.5 or ^ ^ 30° fe.g. lSuenzh and Schmi J2009( : ISchmid et al. 



2006al ) which is roughly the location of the equator in our observations. A slightly higher 
lim b polarization is possible with a stratospheric haze model similar to the one determined 
by iKarkoschka and Tomaskd (120051 ) consisting of a relative thick layer with high albedo, 
highly polarizing, and forward scattering particles 

However, the interpretation of the strong equatorial polarization in Saturn as higher 
order scattering or "limb polarization" effect is in contradiction with the observations. Limb 
polarization would produce a Qr/I signal along the entire equatorial band with maxima 
above 5 % at the equatorial limbs. This is not observed. In addition we would expect a 
high polarization of at least p > 20 % at large phase angle s whic h is also not visible in the 
contemporaneous Cassini polarization image of I West et al.l (120091 ). 



An alternative explanation could be strong light scattering from the inclined ring system. 
But this can be rejected with an order of magnitude estimate. The solid angle subtended 
by the illuminated ring system as seen by the equatorial region is small and the forward 
scattering of light by the ring is inefficient b e cause the ring is made of large (> 1 mm), back- 
scattering bodies. For example iDones et al.l (Il993l ) show that the phase function for forward 
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scattering by the ring is only I{a) < 0.1/(0) for a > 135°. Therefore the irradiation with 
hght scattered from the rings is well below 1 % of the direct irradiation from the sun and it 
is impossible that this low level of indirect light produces a significant fractional polarization 
at the equator. 

Which other effect could explain the transient equatorial polarization of Saturn if phase 
effects, limb polarization, and ring reflection can be discarded? The unknown effect seems 
to be restricted to the equatorial band and photon incidence and emergence angles of ~ 25° 
for the incoming and back-scattered radiation. Such a localized feature could perhaps be 
formed by a higher order scattering effect in an optically thick, but geometrically thin cloud 
or aerosol layer. We are not aware of investigations of such scattering geometrie s. However, 



we m ay hope that a detailed analysis of the Cassini polarimetry presented in I West et al. 



(120091 ) clarifies the case. 



6. Summary, discussion, and conclusions 

This work presents ground-based spectropolarimetry and imaging polarimetry of Jupiter 
and Saturn. The new data are of unprecedented quality because modern instruments with 
CCD detectors were used. This type of data opens up many new avenues of polarimet- 
ric investigation for planets, because the polarization signal can be quantified much more 
accurately from imaging polarimetry and spectropolarimetry, when compared to previous 
aperture polarimetry. Well calibrated data with quantified seeing effects allow for detailed 
comparisons with model calculations of the scattering layers in Jupiter and Saturn. 

There are a series of important results in this work: 

We present for the first time spectropolarimetry for the strongly polarized poles of 
Jupiter. The polarization shows an overall decrease with wavelength and a rich spectral 
structure with enhanced fractional polarization in strong methane absorptions. In the po- 
larization flux spectrum the absorption bands are still visible as absorptions, although with 
a lower equivalent width than in the intensity spectrum. The enhancement in fractional 
polarization Q/I m the methane bands can be explained by a highly polarizing stratospheric 
haze layer above an atmosphere where multiple scatterings (or the reflection of unpolarized 
light) are suppressed in the methane absorptions. 

The fractional polarization at lower latitudes of Jupiter is small and slightly negative 
with an enhanced polarization signal in strong methane absorption. 

From the polarization profiles along the central meridian of Jupiter we can derive a 
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latitude of ±58° for the transition between the polar region with strong positive polarization 
and lower latitudes with a negative polarization. Our results agree very well with the polar- 
ization measurements from the Pioneer missions (jSmith and Tomaskdll984[ ). Stratospheric 
haze is responsible for the strong polar polarization and the sharp transition is pointing to 
a well defined border in the stratospheric circulation. 

For the poles of Jupiter we measure in the V-band a seeing limited peak limb polarization 
of +9 to +10 %. This indicates a resolution-corrected peak polarization of about +11.5 %. 
All previous studies reported lower values (p ~ 6 — 8 %) for the maximum polarization at 
the poles of Jupiter. This difference can be explained by the better spatial resolution (seeing 
~ 1") of our data and the appropriate seeing correction. 

We derive the polarization flux of the entire positively polarized polar hoods (called S+ 
and N+ in this work). This approach should allow an easy comparison of our data with 
previous and future observations for long term studies which are interesting for investiga- 
tions of the haze production, destruction and transport in t he polar stratosphere of Jup iter. 
Such long terrn polar ization changes were reported e.g. by IStarodubtseva et al.l (120021 ) and 
Starodubtseval feoool ). 



We found no limb polarization models with a limb polarization higher than 10 % in 
the literature. For this reason we carried out exploratory calculations. They indicate that a 
polar haze l ayer consisting of f orwar d scattering and highly polarizing aggregate particles as 
proposed by IWest and SmithI ( 119911 ) is compatible with our observations. 



The polarization of Saturn is lower than that of Jupiter. Our observations of 2002 and 
2003 revealed at 6000 A a positive feature (about 1.4 %) at the South pole and a stronger 
(> 1.6 %) transient polarization feature near the equator. 

Spectropolarimetrically, Saturn shows hardly any enhancements in the fractional polar- 
ization at wavelengths of strong methane bands in strong contrast to the case of Jupiter. 
This indicates that the polarizing particles in Saturn are not located in a discrete layer, but 
well distributed in the scattering atmosphere. Absorptions within the polarizing scatter- 
ing layer reduce the reflected intensity / and polarization flux Q simultaneously and causes 
therefore no or only weak features in the fractional polarization spectrum Q/I. 

At the South pole and the equator the polarization decreases rapidly with wavelength, 
and essentially no polarization p < 0.4 % is present for wavelengths above 7500 A. The fast 
decrease of the limb polariz ation with wavelength is indicati ve for small scattering particles 
(< 100 nm) as proposed by iKarkoschka and Tomaskd (120051 ) based on HST imaging. 



The strong equatorial polarization feature in Saturn is remarkable since it was not 



-37- 



present in earlier measurements. During our observations the inclination of Saturn was high 
and therefore the effect could be seasonal. We investigated various possibilities to explain the 
equatorial polarization but could not find a reasonable solution. Despite this we are certain 
that the observed fractional polarization of about p = 1.6 % in polar direction is real. It 
seems interesting to investigate this feature further with other data or model simulations. 

An important conclusion from this paper is that modern polarimetric measurements of 
Jupiter and Saturn from the ground can provide accurate quantitative results which constrain 
strongly the scattering properties of the atmospheres. This type of investigations certainly 
did not receive sufficient attention in the last decades, and more observations, e.g. for other 
wavelengths, and a lot of limb p olarization modelling s till ne eds to be done. A first step with 



model simulations was done by iBuenzli and SchmidI (120091 ) who present an extensive grid 



of model calculations for Rayleigh scattering atmospheres. Further models ought now to be 
calculated which should investigate the dependence of the limb polarization signal on the 
scattering phase matrix of different populations of haze particles. In additions one needs to 
take into account the spectropolarimetric structure in the methane bands in order to derive 
the vertical stratification of the scattering layers. 

Scattering layers, refiecting the solar light, are an intriguing aspect of solar system 
planets. They affect the radiative transfer in these objects. For the investigation of the 
reflected light from extra-solar planets a comprehensive understanding of the physics of the 
high altitude haze layers is very important. For this reason it is essential to carry out 
detailed investigations of the reflecting layers in solar system planets. Investigations based 
on modern polarimetric observations, as presented in this work for Jupiter and Saturn, are 
therefore very valuable for progress in this direction. 
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